COMPLETELY REGULAR MULTIVARIATE STATIONARY 
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Abstract. We are going to give necessary and sufficient conditions for a multivari- 
ate stationary stochastic process to be completely regular. We also give the answer 
to a question of V.V. Peller concerning the spectral measure characterization of 
such processes. 



1. Introduction 

In this paper we shall give a necessary and sufficient condition for a multivariate 
stationary stochastic process to be completely regular. For the scalar case the descrip- 
tion of completely regular processes was obtained by Helson an Sarason, see @, ^. 
Almost none of the scalar methods is available in the vector situation. The explana- 
tion is simple. Our problem will be reduced to verifying weighted inequalities for 
a certain integral operator. The weight will be a matrix weight arising from the spec- 
tral measure of the process. All the pointwise estimates of integral operators become 
too crude for the vector valued case. For example, if a positive kernel is majorized 
by another one, and this second kernel gives the bounded operator in L'^{fi), then the 
original kernel obviously corresponds to a bounded operator in L^{fi) too. But this 
is not the case if /z is a matrix measure even for scalar kernels. 

The study of prediction theory for multivariate stationary stochastic processes was 
started by Kolmogorov and Wiener in the 50's, see, for example |T^, [Q, and 0]. It 



was later continued in works of I. Ibragimov, Yu. Rozanov, V. Solev, A. Yaglom, V. 
Peller, S. Khruschev, N.J. Young. An extensive bibliography can be found in |^ (for 
scalar processes) and in (for vector ones). 153al58,169 

Let us recall that a multivariate stationary stochastic process with discrete time 
is a sequence of (i-tuples x{n) = {xi{n),X2{n), . . . ,Xd{n)), n G Z of scalar random 
variables such that E|a;j(r;,)p < oo and the correlation matrix Q{n, k) 



Q{n, k) = {Q{n, A;)ij}^<. := <^ Exi{n)xj{k) 

— ■J — \. ) l<i,j<d 

depends only on the difference n — k; here E denotes mathematical expectation. 

It is well known (see [§]) that there exists a matrix- valued non- negative measure M 
on the unit circle T whose Fourier coefficients coincide with entries of the correlation 
matrix 

Q{n,k) =M{n-k), n,keZ. 
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The measure M is called the spectral measure of the process {x{n)}n<^x- 

The random variables Xj{n) can be treated as elements of Hilbert space L'^{Q, dP), 
where Q is the probability space and P is the probability, so x{n) can be treated 
as elements of the M'^- valued space L^^(f2,(iP) For a moment n of time we can 
consider the past Xn and the future A"" of the process, which are defined as the 
subspaces 

Xn = span {xj{k) : 1 < j < d, k < n} 
A"" = span {xj{k) : I < j < d, k > n} 

of L\Q,dP). 

A process is called regular if n„>oA:'" = {0}. In this case (see 0) the spectral 
measure M of the process is absolutely continuous with respect to Lebesgue measure. 
Let W be the density of M with respect to Lebesgue measure. The matrix-valued 
function W is called the spectral density of the process. 

A process {x{n)}n<^z is called completely regular if its past is asymptotically or- 
thogonal to the future, namely if 

sup {|E(^r7)| : ^eXo,r]^ < 1, E|r/|2 < l} — ^ as n ^ oo . 

Of course, complete regularity implies regularity. If the process is Gaussian (i.e. all 
random variables Xj{k) have normal distribution) then the complete regularity means 
simply that past and future are almost independent. The problem we are dealing 
with is to characterize completely regular processes in terms of spectral measure. 

It is well known (see again ||^) that if the process is completely regular, then its 
spectral measure is absolutely continuous, dM = Wdm where dm is the normalized 
(m(T) = 1) Lebesgue measure on the unit circle T. 

The reader is referred to once more to see that there exists d^ < d (the rank of 
the process) such that the spectral density W{t) has rank d^ for almost all t G T. If 
do = d then the process {x{n)} is said to be a full rank. 

The study of processes of arbitrary rank can be easily reduced to the study of the 
processes of full rank, see [Q . So in this paper we shall consider only processes of full 
rank. 

For the scalar case the description of completely regular processes was obtained by 
Helson an Sarason, see @, |^ . To state their result we need a couple of definitions. 

Let us recall that a function / on the unit circle T belongs to the space BMO 
(bounded mean oscillation) if 

sup ^ 1/ - fi\dm = II/IIbmo < oo ; 

here // denotes the mean value of / on the interval /: // := Jj fdm and the 
supremum is taken over all subarcs I of T. 

The space VMO (vanishing mean oscillation) consists of all function / G BMO 
such that 
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Theorem 1.1 (Helson, Sarason). Let w be the spectral density of a scalar sta- 
tionary process. Then the process is completely regular if and only if w admits a 
representation 

w = \p\'^e'^ , 

where p is a polynomial with roots on the unit circle T and (f is a real-valued function 
in VMO. 

It was conjectured by V. Peller in P| that the same resuh holds for muhivariate 
stationary processes. Namely he conjectured that a multivariate stationary process 
is completely regular if and only if its spectral density W admits the following rep- 
resentation 

W = P*e'^P, 

where P is a polynomial matrix whose determinant has roots on T and the matrix 
function $ = $* belongs VMO. 

In this direction he was able to prove the following theorem 

Theorem 1.2. A multivariate stationary process is completely regular if and only 
if its spectral density W admits the factorization 

W = P*WiP, 

where P is a polynomial matrix whose determinant has roots on T and Wi is the 
density of a completely regular stationary process such that E L^. 

1.1. The main result. Let us recall that a measure fi on the unit disk D is called 
Carleson if 

sup /i(g(/)) <c-\i\ 
I 

and is called the vanishing Carleson measure if 

limsup/i(Q(/))/|/| = 

|7H0 

where limsup is taken over all subarcs / of T. Here Q{I) denotes the "Carleson 
square" for the arc I, 

Q(I) = {ze'B : z/\z\ El,l-\I\< \z\ < 1} 

For a function F on the unit circle let F{X), A G D, denote its harmonic extension 
at the point A. 

The main result of the paper is the following theorem. 

Theorem 1.3. Let the density W of a stationary process satisfy e . Then 
the the following are equivalent 

1. The process is completely regular; 

2. is the spectral density of a completely regular process; 



3. limsup 

all subarcs I of T; 



/ I r \ 1/2 / 1 /■ \ 1/2 



1; here supremum is taken over 
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1/ 



W{\)) (W'\X) 



4. limsup 

|A|-1 

extensions of functions W 



1/2 



1, where W{X) and W ^{X) are harmonic 



T and W' 



5. limsup|det(^H^(A)j exp(^- [logdet VT] (A) 



are 



harmonic extensions of functions W 

AG©. 

6. The measures 



T respectively at point AG©. 

1, where W{\) and [logdet W] (A) 

T respectively at point 



T and log det W 



d 



W{z)-^/^ ^W{z) W{z)-^/^ 



and 



W{z) 



-1/2 



dx 

d_ 

dy 



W{z) W{z 



-1/2 



\z\^)dxdy 



'1 — \z\^)dxdy 



are vanishing Carleson measures 



Together with Theorem |1.2| the above theorem yields the complete description of 
completely regular stationary processes 

Theorem 1.4. A stationary process with spectral density W is completely regular 
if and only if W admits the representation 

W = P*WP, 

where P is a polynomial matrix whose determinant has roots on T and the matrix- 
function W satisfies W^^ G and one of equivalent conditions 3-6 of Theorem 

M 



Let us discuss the main result (Theorem |1.3| ) a little bit. First of all it is not 
difficult to show directly that in the scalar case the conditions 3-6 of Theorem [L3 
are equivalent to = e*^, G VMO. We are leaving this as an exercise for the 
reader. 

Usually in probability only real valued stationary processes are considered. In that 
case the spectral density of a process should satisfy W{j) = W{z), and only such 
functions can be realized as densities of stationary processes. 

If one allow complex-valued processes, any non-negative matrix function is the 
spectral density of some stationary process. 

Our theorem deals with arbitrary non-negative matrix-functions and can be applied 
to complex- valued processes (as well as to real- valued) . 



2. Scheme of the proof of the main result 

The diagram of the proof will be the following: 1 4 5 => 6 
we will show that 1 2 and so automatically 2 =^ 1. 
And in this section we will show that 3 <^==^ 4. 



1. Then 



Lemma 2.1. For a scalar weight w the following conditions are equivalent: 
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1. limsupf TjT [ w] (— I w ^ ) = 1; 



2. limsup w(A)w~"'^(A) = I, where w{\) and w~^{\) denote the harmonic exten- 

|AH1 

sions of w and respectively at the point A; 

3. w = e^, where G VMO. 

Proof. First of all let us rewrite condition 1. Let if := logw. For a function / let 
denote its average over the arc I, f^ := |/|~^ /• Then clearly 

^/ ■ = h/exp(-y?^)] ■ [(w-^)^exp(¥?^)]. 

By Jensen inequality (geometric mean < arithmetic mean) the expressions in brackets 
are at least 1, so the condition ! splits into the following 2 conditions 

lim sup [w, exp =1, and limsup [w""*^ exp =1. 
|/Ho |/Ho 

Let /+ denote the positive part of the function /, := max(/(x),0). Then the 

inequality 

X <e^ — I for X > 

implies 

\l\ " "^^^^ ~ Wlli^^^^^'^ ~ ~ = ^/ exp(-^^) -1^0 

Since fj\ip — ip^l = 2 fj{f — fj)+, one can conclude that ip G VMO. 

Similarly, using Poisson averages instead of averages over intervals one can get 
from condition 2 of the lemma that harmonic extension of \(f — <f{X)\ at the point A 
tends to as A ^ 1. But that is an equivalent definition of VMO, so the condition 
2 also implies that ip G VMO. 

On the other hand, if G VMO, John-Nirenberg Theorem (see |l| [Chapter VI], 
the measure of the set {t G / : \'p(t) — 'Pj\ > a} is estimated from above by Ce~^", 
where K = Kj — » oo as |/| ^ 0. Therefore for x > 1 the measure of the set {t G / : 
exp{(p(t) — ipj) > x} is estimated from above by Cx^^ . Integrating this distribution 
function one can get that limsupiJi^o""^/ 6xp(— < 1 (in fact, it is 1, because 
by Jensen inequality w^exp(— (/j^) > 1. Similarly, limsup|j|^g(w~^)^ exp((y9j) = 1. 
Multiplying the above two inequalities one gets condition 1. 

The proof that 3 =^ 2 is similar. For a point A G D let Ix be an interval with center 
at A/|A| of length — Since the Poisson Kernel Px{z) = (1 — |Ap) ■ |1 — \z\'''^ 
satisfies sup^g^w Pa(-2) ^ as |A| 1, the distribution inequality for ip on / 



as 1/1 — > 0. 



A 



implies that w{X) ■ exp(— y9(A)) — > as |A| 1, and therefore the condition 2 of the 
lemma. □ 

The following Lemma is probably well known and can be easily from the distribu- 
tion function inequality for VMO (John-Nirenberg Theorem). 

Lemma 2.2. For X & 3 let Ix be an interval centered at A/|A| of length 1 — |A|. 
Ifpe VMO, then p^ - p{\) as |A| 1. 
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Corollary 2.3. Let (p G VMO and let w = e'^ . Then for Ix as in the above lemma 
we have 

ii,„ ^ = 1. 



|A|^l W 



Proof. By the above lemma lim|A|^i exp(w(A))/ exp{Wj ) = 1. On the other hand it 
follows from the proof of Lemma O] that 



hm w{\)/ exp{ip{X)) = 1 and hm w^/ exp{(pj) = 1. 

Taking the ration of the last 2 identities (with / = I\)we get the statement we 
need. □ 



Now to show equivalence of condition 3 and 4 of the Theorem |L3| is enough to 
show that these conditions imply that for a fixed vector e G scalar weight w{z) = 



(W{z)e, e) satisfies conditions 1 and 2 of Lemma Then Corollary |2.3| implies 
that the averages Wj^ and W{X) are equivalent, the same holds for W~^, and we are 
done. 

It remains now to show that the scalar weight w{z) = (W{z)e, e) satisfies condition 
1 (equivalently 2) of Lemma |2.1| . The easiest way to do that is to recall where the 
Muckenhoupt condition {A2) came from, see ||10|| . 

Recall that the quantity || is just the norm of the averaging 
operator / ■ in the weighted space L'^iW), see [p!0| [Lemma 2.1]. Then 

[tt;^]^/^[(w"^)j-]^/^ is the norm of the restriction of the above averaging operator onto 
the subspace of L'^{W) consisting of functions of form fe where / is a scalar function. 
Therefore 

so the weight w satisfies condition 1 of the lemma. 

Similarly, the quantity ||iy(A)^/W-^(A)^/2|| jg j^^g^ ^-^^ 

norm of another averaging 

operator (/ 1-^ Jr^fkxjkx, where kx is the normalized reproducing kernel of H^, 
kx{z) = (1 — |Ap)^/^(l — A^;)""*^, see [Lemma 2.1], so condition 4 of the theorem 
imply condition 2 of the lemma for the weight w. 



3. Eliminating probability 

The problem of description of completely regular processes can be now stated 
without mentioning any probability theory at all. 

First of all notice that without loss of generality we can assume that the process 
is complex- valued. Namely, if we have a real stationary process {x{n)}n(^z we can 
consider its comlexification, namely the same process but in the complex Hilbert 
space L'^a{^,dP). Considers the comlexificated past ('^n)^ future i^")^ 

{'^n)c = span {xj{k) : 1 < j < d, k < n} 
(A:'")c = span {xj{k) : 1 < j < d, k > n} 
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where span now means the closed linear span in the complex Hilbert space L'^^i^, dP). 
It is easy to see that 

sup {|E(er7)| : ^ G A-o, G X'\ E|eP < 1, Elr^l^ < l} = 

= sup {|E(er7)| : i G (A'o)c, V e ('^")c, < 1, < l} , 

so a process and its comlexification are completely regular simultaneously. So we 
indeed can assume from the beginning that our process is complex valued. 

Consider now the vector space L'^{W) of C^- valued functions on the unit circle 
with the norm 

(of course we have to take the quotient space over the functions of norm 0). The 
mapping Xj{k) i— > z'^ej, where ej, j = 1, d is the standard orthonormal basis in C^, 
is an isometric isomorphism between spaja{xj{k) : 1 < j < d, k & Z} and L'^{W). 
The past Xn and future X"' are mapped to the spaces X„ and X" of L^{W) 

Xn = spanj^'^C'^ : k<n} (3.1) 
X" = spanjz'^C'^ : k > n} . (3.2) 

So the problem of describing completely regular stationary processes can be refor- 
mulated as follows: describe all matrix weights W such that the spaces Xq and X" 
are asymptotically (as n ^ oo) orthogonal to each other, 

Pn = sup{|(e,r/),,(^^| : e e Xo, G X^ H^H^.^^^ < 1, \\v\\,.^^^ < l} ^ 0, 

(3.3) 

as n ^ oo. 



4. Necessity (1 =^ 4) 



In this section we are going to prove the implication 1 
below) and the equivalence 1 <^==^ 2 (see Lemma 



4 (see Theorem 4.1 



For a function F defined on the unit circle T let F{X) denote its harmonic extension 
a the point A G D. 

Theorem 4.1. Let W be a matrix valued weight such that G . Suppose 
the "past" Xq and "future" X" defined by (|3.1| ), (|3.2|) are asymptotically orthogonal, 
that means 



pn 



as n 



supi |(^,r/) 



l2(vi/)| 



: e e Xo, e X", lieil,,^^^ < 1, < 1 







oo. Then 



lim sup 

|A|-1 



w{\)y\w-\x)) 



1/2 



1. 
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Proof. First of all let us show that if ^ is completely regular and W ^ E then 
W satisfies the Muckenhoupt (^2) condition 



sup 

A6D 



w{x)y^^(w-\x) 



1/2 



< 00 . 



Recall that 

{f,kx)kx in t 
ing kernel for H^, 



wi\)y^^(w-\x) 



1/2 



is exactly the norm of the operator / 



le weighted space L'^{W); here kx denotes the normalized reproduc- 



kx{z) := 



|A|2)l/2 



l-\z 

\\kx\\2 = 1- Note that ko = 1. So if W^^ G the operator / 1-^ (/, 1)1 is bounded in 
L'^{W), and therefore by translation invariance the operators / 1-^ (/, z^)z'^ = f{n)z^ 
are bounded as well (they all have the same norm). 

We know that the spaces Xq and X"- are asymptotically orthogonal, so we can say 
that for large enough N the operator P+ restricted onto span{Xo, X'^} = span{;2"'C'^ : 
n ^ [0, N]} is bounded, say by 2, 

V/ e span{Xo, X^} = span{;s"C'^ : n ^ [0, N]} . 



Since / - En=o/(^)'^'' ^ span{Xo,X^} = span{z"C'' : n ^ [0,N]}, one can con- 
clude that the operator P+ is bounded in L'^{W) and so the weight satisfy the Muck- 
enhoupt condition (A2). 

We will need the following simple lemma about Muckenhoupt weights. 

Lemma 4.2. If w is a scalar Muckenhoupt weight, then its harmonic extension 
w{X) can't decay to fast near the boundary of the disk. Namely, if the Muckenhoupt 
norm of w is at most M there is a function a — a^^, a : [0, 1) (0, 00), a{t) \ 
as t ^ such that 

(1- iAnw(o) 



w{X) 



<a{\X\). 



Proof of the lemma. For an arc 7 C T and A; > let kl denote the arc of length k\I\ 
with the same center as /. 

We are going to show that for a Muckenhoupt weight w with the Muckenhoupt 
norm at most M 



,j <M\2- e)''wj , e = e{M) > . 



(4.1; 



Applying this formula in the case 2"/ = T and using the trivial estimate 

w(A) > Cwj^ 

where Ix is the arc with center at the point A/|A|, |Ja| = 1 — |Ap and C is an absolute 
constant, we can get from there (recall that — 1 — |Ap = 2~") 

w{X) > c(2 - £)-" • w{0) = c(2 - £)i°g2(i-|Ap) . ^(0) = c • (e - 5)'°g(i-l^l') . w{0), 
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where 6 = 6{6) > 0; here e is the base of the natural logarithm, not a vector in 
C^. This estimate implies the conclusion of the lemma with a{t) = c^^(l — t^) ■ (e — 

To prove ( |4.1| ) we notice the since the weight w"^ is the Muckenhoupt {A2) weight 
with the same Muckenhoupt norm as w, it is doubling and therefore 

where e depends only on the Muckenhoupt norm of w. Iterating this inequality n 
times we get 

The last estimate and the Muckenhoupt condition imply 

< M/{w-')^„j < M ■ (2 - er/{w-')j < M^Wj 
and that is exactly what we need. □ 

Corollary 4.3. // a matrix weight W satisfies the Muckenhoupt condition {A2) 
with the Muckenhoupt norm at most M then for any e G C"' 

where a = um is the function from Lemma \4-2^ - 

Proof of the corollary. The proof follows immediately from the fact that the scalar 
weight w, w{^) = (W{^)e,e)^ is the Muckenhoupt {A2) weight with the Mucken- 
houpt norm at most M (see fill], proof of Corollary 2.4). □ 

We now return to the proof of the theorem. 

The condition e implies that J.^ log detW{^) dm > —00, hence (see [0) 
there exists a factorization of W of the form W = F*F, where F is an outer matrix 
function in H^. 

Take e G and let us compute the distance 

^^^*L2(^^){^"^e,span{;z"C'^ : n > 0} = dist^2(t^){e, spanjz^C^ : n > 0} . 

By the vectorial version of the Szego theorem (see 0) this dist ance is exactly ||-F(0)e||. 
Using the Mobius transformation of the disk one can get from there 

^^^HW)^ J- X ^^^P^^{^"^' : n > 0}} = ||F(A)e||^, . 
Writing the Fourier series expansion of — — 



|\|2Nl/2 °o 

\2U = (i- lA^i/^ J2 A"^-("+^) 



z 

n=0 
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one can see that for any fixed N > the function ^l^^— e is almost in the "past" 
as |A| 1. Namely, 

z — A 

n=0 n=N 

where /| G and /| is small, 

as |A| — * 1, where a{.) is as in Lemma [4.2| and Corollary [4.3| . 

Since Xq and are asymptotically orthogonal, the shift invariance implies that 
the subspaces X_7v and X*^ are asymptotically orthogonal as well. Taking |A| — 1 
and then X — » oo we can conclude that 

\\nX)e\\^J\\W{XY/'e\\^, = 
= d\st^,^^ /^~Jf^'^\ ,span{z-C' : n > 0}}/ ||M|^^^^^ > 1 - /3(|A|)V2 - 1 , 

where depends only on the Muckenhoupt norm of W and /?(|A|)— >Oas|A|-^l. 
The last inequality implies 

\\W{Xy/'F{X)-'\\<{l-P{\X\))-K (4.2) 
Note that since \\F{X)e\\^J\\W{Xy/^e\\^^ < 1 for all e G C^, we have 

\\W{Xy/^F{X)-^\\ > 1. 

We will show a little later that under assumptions of the theorem the subspaces Xq 
and X^ in the weighted space L'^{W~^) are asymptotically orthogonal as well. The 
factorization W = F*F yields the factorization = of W~^. Similarly 

to the previous case 

dist^,^^_,^{^i^l^£^e,span{z-C'^ : n > 0}} = ||F-i(A)*e||^, = \\FiX)-'*e\\^,. 
Acting as before we get 

\\W-\Xy/'F{Xy\\<il-(3,{\X\)r' (4.3) 

where /5i(|A|) — > as |A| 1. 

Combining ( ^7^ ) and ( [4. 3] ) we get 

\\w{xy/'w-\xy/'\\ < (1 - m\)r\i - mM)r' - 1 as |ai ^ i . 

So, we completed the proof modulo the following lemma. □ 
This lemma also gives us the equivalence 1 <^=^ 2. 
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Lemma 4.4. Under assumptions of Theorem \4-.l\ the weight is a spectral den- 
sity of a completely regular process, i. e that the spaces Xq and are asymptotically 
orthogonal (as N ^ oo) in the weighted space L'^{W~^). 

Proof. It is enough to show that 

||P+ I span{Xo,X^}||^,^^_,^^^,^^_,^ ^1 as iV ^ oo. 

The later is true because 

||P+ I span{Xo,X^}|U_,. = \\W-'/'{P+ \ span{Xo, X^})!^!/^,, 



= \\W'/'{P^ I span{Xo,X^})W^-V2||^^^^^ = ||p^ | span{Xo, X^}||^,^^^^^,^^^ 
and 

||P+ I span{Xo,X^}||^,^^_,^^^,^^^ ^1 as iV ^ oo 
(since Xq and X^ are asymptotically orthogonal in L'^[W)). □ 

5. Vanishing Carleson measures 

Recall that Vr(A) and W^^{X) denote harmonic extensions at the point A G D of 
the weights W and respectively. 

Lemma 5.1. Let a matrix weight W satisfy 



hm\\Wi\Y/'{W){\y/'\\ 



Then 



limsup|det(^Vr(A))exp(^-[logdetl^](A))} = 1. 



Proof. First of all let us notice that the assumption of the lemma implies that 
W,W~^ G L^(T), therefore log(det ly) G L^(T). Therefore there exists a factor- 
ization W = F*F a.e. on T, where F is an outer function in H'^i^M^ixd) 

Since F is an outer function in iY^, det F is an outer function in H'^^'^. Therefore 

|detF(;z)| = exp{(log|detF|) (z)} = exp | i (log det 14^) (2) | (5.1) 

It is well known fact that F*{z)F{z) < W{z) for any z ^ D, where < means the 
inequality for quadratic forms. There are many proofs of this fact, for example it 
admits a very simple operator-theoretic interpretation which is in fact hidden in the 
proof of Theorem Explanation that we present here is more function-theoretic: 
Direct computation shows that 

A {F{zyF{z)) = 4 {dF{zy) {dF{z)) = 4 idF{z)y {dF{z)) > , 

so for any e G the function ||F(2)e|p is subharmonic and coincide with (W{C,)e, e) 
on T. 

We can do the same factorization for W~^. Namely, let G be an outer matrix- 
valued function in H'^{Mdxd) such that = G*G on T. We should point out 
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to the reader that in general G does not necessarily coincide with F^^. However, 



applying ( |5.1| ) to G one can conclude that 

|detG'(;z)| =exp|^ (logdetiy^^) = \detF{z)\-^ 

Now we are in position to prove the lemma. By the assumption 

1, 



\imJW{zf'\W-^){zY/^\ 



(5.2) 
(5.3) 



and therefore, 



lim^ \^ei{W{z)) det {{W-^){z)) \ 



Using ( |5.2| ) one can rewrite the last identity as 

lim {[<letW{z)/\detF{z)\^] [<leiW-\z) /\deiG{z)\^]] = 1 



2^1 



Since F{z)*F{z) < W{z) and G{zYG{z) < W ^{z), expressions in brackets are at 
least 1, so, taking into account ( p.l| ) 

lim [det Vr(2)/exp{(logdet Vr)(2)}] =0 



or equivalently 



lim log{det(iy(z))} - (log det PF) (z) = . 

\z\^l 



(5.4) 
□ 



Theorem 5.2. A matrix weight W satisfy 



lim sup I det \W{\)j exp [log det W] (A) 



if and only if the measures 



d 



and 



W{z)-^'^ ( ^W{z) ] W{zY^'^ 



W{z)-^'^ ( ^W{z) ) W{z)-^l^ 



[1 — \z\^)dxdy 



[1 — \z\^)dxdy 



are vanishing Carleson measures. 



The implication 3 =^ 4 of Theorem follows immediately from Theorem |5]2 
and Lemma |5.1| . 

To prove the theorem we need the following well known description of vanishing 
Carleson measures 

Lemma 5.3. A measure fi in the unit disk ID) is a vanishing Carleson measure if 
and only if 

1-IAP . , , 
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We also need the following lemma that was proved in |]rT|, see Lemma 3.1 there. 

Lemma 5.4. Let W be a harmonic function ofn variables with values in the space 
of strictly positive dx d matrices (W{x) = W{x)* > Vxj. Then 

A(log(detW^)) = -^trace ({W-^'^^W'^'^f 
j=i ^ 

Proof of Theorem |5.^ . The proof below follows the lines of the proof of Theorem 3.2 
of 0. 

By Green's formula and Lemma 



1 

2^ 



// log 


1 — sz 


J J o 


z — s 



^// (trace 



.mdWiz) 



A log {det {W{z))} dxdy 

2 



dx 



W{z) 



A/2 



+ 



-mdWiz) 



+ trace ( W{z)-'/^^^^;^W{z)-'/^ 

oy 



lot 



1 — sz 



z — s 



dxdy 



Using an elementary inequality log(l/a) > 1 — a for < a < 1 and the fact that 
< traced for a non- negative matrix A, the last integral is at least 



1 

47r 



W{z) 



dx 



1 



Together with ( |5.4| ) this imply 



W{z^-'^ 



log 



1 — sz 



Wizr^^^Wiz)-'/^ 



z — s 

2 



dxdy > 



dx 



-mdWiz) 



dx 



W{z 



1 



1-1/2 



1 — sz 




j dxdy = 


z — s 


' (1- 




1 — sz\^ 



dxdy 



lim 



(1- 


s\ 




ll -■ 


sz 


2 



w{z)-y-^-^w{z)-y- 



dx 



(1 — \z\^) dxdy = 



that yields that the measure \\w{z)-'/^ {£w{z)) w{z)-'/^\\^ (1 - i^n dxdy is a 

vanishing Carleson measure. 

W{z)-^/^ W{z)-^/^ (1 - dxdy is treated similarly. 



The measure 



To prove the opposite implication, let us estimate the integral 



trace ( W{z)-'/^^^^W{z)-'/' ) log 



1 — sz 



z — s 



dxdy 



(the integral with dW/dy can be estimated absolutely the same way). Denote by bs 
a Blaschke factor with zero at the point s, bs{z) = (z — — sz)^^. 
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First of all, we can estimate the trace by d ■ || ■ ||, where d is dimension of the space. 
So we can estimate the integral by 



c 



W{z) 



dx 



\og%{z)\-^dxdy= JJ ■■ + J J 



\bs{z)\<S 



\bs{z)\>e 



To estimate the second integral we notice that 
\og\bsiz)\-^dxdy < C{e 



(1- 


s 






z 






1 


— sz 


2 



for > e, and since the measure is a vanishing Carleson measure we can make 

the integral as small as we want when \s\ 1. 

To estimate the first integral let make a trivial observation: if w G L^{T), w > 
and w{z) denotes its harmonic extension at the point z, then for all z such that 
\z\ < 1/2 (and therefore for all z such that 1^1 < e < 1/2) 

^wiz) < CwiO) , 

where C is an absolute constant. Combining this observation with the Harnack 
inequality w{0) < C'w{z), \z\ < 1/2, and applying it to functions w{.) = (W{-)e, e)cd 
we get the inequality for quadratic forms 

d 

— < CeW{0) < CiW{z) . 
ox 



It in turn implies 

^^'^(^W{z)]W{z)-^'^ <Ci, \Jz: \z\<e<ll2. 
\dx / 

Using the Mobius transformation z i— > hs{z) we get 

d 



W{z)-^'^(^-^W{z)^W{z)-^/^ 



<Ce, yz: \bs{z)\ <e<l/2. 



Since 



1 



log \ bs{z)\ dxdy < Ce log - , 

\bs{z)\<e ^ 

we can estimate the first integral by Ce^ log(l/£); we can make this number as small 
as we want by picking sufficiently small e. 

□ 



6. Embedding theorem and equivalent norms 



By analogy with the scalar case (see [jT2[) we will say that a matrix weight W 
satisfies the invariant Ar^ condition if 



sup <^ det ( W{s) ) exp ( - [log det W] {s 



< oo . 



{invAc 



The supremum is called the invariant A^o norm of W. 
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Theorem |5.2| implies that if the measures 

2 

|2\ 
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-1/2 (^^W{z)^ W{z)-^/^ 



and 



d 



w{zy^'^ ( -^w{z) ] w{zY^'^ 



(1 — \z\ )dxdy 



\z\'^)dxdy 



are vanishing Carleson measures then the weight W satisfies the invariant A^o con- 
dition. 

Literally repeating the proof of Theorem |5.2| one can obtain that the weight W 
satisfies the invariant A^o condition if and only if the above measures are Carleson. 
We will need the following "embedding theorem" . More general result was proved 

in 



11 , Lemma 4.1. 



Lemma 6.1. Let W be a matrix weight satisfying the invariant Aoo condition, and 
let II he a Carleson measure with the Carleson norm ||yu||^. Then for any analytic 
(or antianalytic) vector-function f , the following inequality holds, 

{W{z)f{z), f{z)) df,{z) <CM\J {W{i)f{i), f{0)dm{0 , 

where the constant C depends the dimension d and the invariant A^o norm ofW. 

Proof. The invariant A^o condition implies that logdetVT G L^, so there exists (see 
0) an outer function F G H'^{Mdxd) such that W = F*F. It is well known (see 
again M) that 



det F{z)\= exp | ^ [log det W] (z) | . 



It is well known and it was already shown it in the proof of Lemma xl that F(z)*F{z) < 
W{z). Hence 

\\W{zy^^F{zy^e\\>\\el e G 

Since 

1/2 



(6.5) 



\det {W{zy/^F{z)-^}\ = |det(^W^(A))exp(^-[logdetVr](A) 



< C 



we can estimate 

\\W{zy/^F{z)-'e\\ <C. 

Together with ( |6.5|) it implies that (W{z)e,e) and ||F(z)e|p are equivalent in a sense 
of two-sided estimate. Therefore 



{W{z)f{z), f{z)) d^{z) <C {F{z)f{z), F{z)f{z)) dfx{z) < 



< CM^. / (F(0/(0, F{Of{0)dm{0 = CM^ {WiOfiO, f{0)dm{0 ■ 



□ 



We also need the following simple lemma. 
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Lemma 6.2 (equivalence of weighted norms). Let W be a matrix weight satisfy- 
ing the invariant A^o condition. There exist a constant C such that for any analytic 
or antianalytic vector- function f in iF'iW) satisfying /(O) = 

^ ! {Wf, f)dm < [[ iWiz)f'iz), f'{z)) log ^ dxdy <C f (Wf, f)dm 
Proof. Let us recall the the operators d and d are defined as 

5/ = 1 f ^ _ M) df = -(^+t^ 

2 \dx dy J ' 2 \dx dy 

Recall that for analytic functions df = f and Of = 0. 

Let / be an analytic function, /(O) = 0. Using the Green's formula and taking 
into account that /(O) = and A = Add = Add we get 

{Wfj)dm=^ [[ A{W (z) f (z) J (z)) log ^dxdy = 

^TT J J p \Z\ 

-If (dW{z)f'{z), f{z)) log ^dxdy +- [[ {dWiz)fiz), fiz)) log ^dxdy 

+ -11 {W{z)nz)J'{z))\og^^dxdy=-{l,+l2 + -I,) 

The last integral X3 is exactly the integral we want to estimate. Let us denote 
:= Jj{Wf, f)dm, B'^ := X3. We want to show that A x i? in a sense of two sided 
estimate. Let us estimate Xi: 



+ 



< 



{W{z)-^/^dW{z)W{z)-^/^W{zy/^f'{z), W{zy/^f{z)) log ^dxdy 
) 1^1 

\Wiz)-^/^dW{z)W{z)-'^^ ■ \\W{zy^^f{z)\\ ■ \\W{zy^^f{z)\\ ■ log -^dxdy 



< (I ^J|H^(^)-^/W(^)W-(^)-i/2||'(iy(z)/(z),/(z))^,log^rfxrfi/^ 

1 1' JW{z)f'{z) J' (z))^ Jog ^^dxdy 



z\ 

1/2 

X 



1/2 



The measure ^/'^dW{z)W{z) ^^"^^ log j^^dxdy is Carleson, so by Lemma O 



the first term in the product is estimated by KA {K is a constant). The second term 
is just B so |Xi| < KAB. Similarly jXsl < KAB. So 

A^ = B^ 

where 

|Xi|, IX2I < KAB. 

This immediately implies 

^A<B<CA 

for an appropriate choice of C. □ 
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To prove the implication 6 =^ 1 we need to estimate Jj{Wz"-f,g)d'm, f G 



Using the Green's formula and taking into account that ^'(0) = and A = Add 
Add we get 



/ {Wz'^f,g)dm = ^ [[ A{Wiz)z-f{z),g{z))^,log^dxdy = 
Jt '^'^ J J o 1^1 

= 1 ff {dW {z)d{zy{z)),g{z))^ Jog ^dxdy + 
J J o 1^1 



+ - // {dW{z){z-f{z)),dg{z))^Jog--dxdy = -(Ji + J2) 

TT / / n, \Z\ TT 



The second integral is easy to estimate: 



< 



< 



{W{z)-^/^dW{z)W{z)-'/^W{zy/^{z''f{z)), Wizy/^dg{z))^, log ^dxdy 

D 1^1 

\\W{z)-^^^dW{z)W{z)-^^''\\ ■ WWizy/^z'^fiz))]] ■ \\W{zy/^dg{z)\\ log -^dxdy 



(l^(z)/(z),/(z))^,log^dxrfyj X 

1/2 



1/2 



X 



/ /^(^^^)'^^^^)'^^(^))c<^ log ^^dxdyj 



The last term is equivalent to the norm ||fi'||^2(^y-) (see Lemma 6^), so by Lemma 6A 



IX2I < ||/||^.(^) ■ ||^?||^.(^) 



■ \\Wiz)-^/^dW{z)W{z)-^/^\\\og^dxdy 



1/2 
c 



Since the measure \\W{z)-^/^dW{z)W{z)-'^/'^\\ log j^My is a vanishing Carleson 



measure, the Carleson norm 



as n — > 00. So IX2I ^ as n — > 00. 



\\W{z)-^/^dW{z)W{z)-^/^\\ log ^My 



1/2 
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To estimate Ti we pick r < 1 close to 1 and split the integral into two: Xi 
/ /rB •••"'"/ /i[i)\rD ' ' ' ' ^cting as with X2 we can estimate 



X 



< 



<l jj ^W{z)-y^dW{z)W{z)-''Y-{W{z)g{z),g{z))^Aog^^^^^^ x 

11 ^{W{z)d{z-f{z)),d{z-fiz)))^Aog^^dxdy^ ' , 



where X is either rD or D \ rD. Note that both terms are uniformly bounded. 

We can say even more. If X = r© the second term can be made as small as we 
wish by picking sufficiently large n. 

Let now X = © \ r©. The measure \W{zY^I'^'dW{z)W{z)-^l'^\\og^^dxdy is a 
vanishing Carleson measure, so for r sufficiently close to 1 its restriction onto D \ rD 



has the Carleson norm as small as we want. So by Lemma |6.1| the first term is as 
small as we want if r is sufficiently close to 1. 

8. A COUNTEREXAMPLE TO PELLER'S CONJECTURE. 

In this section we are going to construct a weight W ^ such that G L^, logiy G 
VMO, but the corresponding stationary process is not completely regular (i.e., the 
weight W does not satisfy any of the conditions 1-6 of Theorem |L 

Let 

\^ 0(2;) y \^ sma cos a 

Here 

5(e^*) = l/log(l/|t|), -l/4<t<l/4, 

and 5 is a continuous function bounded away from and 00 on the rest of the circle, 
and 

«(e^*) = (^/|t|)(5(e^*)^/^ -1/4 < t < 1/4, 

and again a is continuous on the rest of the circle. 
Then 

\ / sin^alog(5 sin a cos a log 5 



\ogW = U* ni e f/-l. ,c 2ir 

y (J log J \^ sm a cos a log cos a log 

and this matrix clearly belongs to VMO: log 5 = log log 1/ It | (considered only in a 
neighborhood of 0) is a "typical" unbounded function in VMO, so cos^ a log S G VMO, 
and all other entries of the matrix are continuous. 

Let us now show that the weight W does not even satisfies the Muckenhoupt 
condition (A2). Direct computations show that 

I cos^ a —sin a cos a \ sin^ a sin a cos a \ 

W = { . .2 ] ■ 2 

- sm a cos a sm a I \ sm a cos a cos a I 
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and 

_ f cos^a — sin a COS « \ f sin^ a sin a cos a ^ 

\^ — sin q; cos a sin^ a j y sin a cos a cos^ a J 

If we pick / to be a symmetric arc [e~*^, e*^] (s > is small), then off-diagonal entries 
of Wj and {W~^)j equal 0, and so we can estimate 

w,>c( '""''^^'^ 
^ ~ \ sm"' a{e) J ' 

^ V ^ "^^^ cos^Q;(e) J 

Therefore 

\\[Wj]^^'^[{W-^]^^'^\\>CS{e)-^sma{e)cosa{e)^oo as £^0. 
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